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Abstract: The main objective of this talk is to develop a matrix pencil approach for the 
study of an initial value problem of a class of singular linear matrix differential equations 
whose coefficients are constant matrices. By using matrix pencil theory we study the 
cases of non square matrices and of square matrices with an identically zero matrix pencil. 
Furthermore we will give necessary and sufficient conditions for existence and uniqueness 
of solutions and we will see when the uniqueness of solutions is not valid. Moreover we 
provide a numerical example. 


1 Introduction 

Linear Matrix Differential Equations (LMDEs) are inherent in many physical, engineer¬ 
ing, mechanical, and financial/actuarial models. In finance for instance, we provide the 
well-known input-output Leondief model and its several important extensions, see [1]. De¬ 
scribing electrical circuits and impulsive behavior is another application described in [1], 
A second order LMDE which describes the vibration of a building is another example, see 
[43]. For other applications of continuous-discrete time systems see [32-35, 44-48]. We 
consider 

A n X^ n \t) + A n _iX^ n P(t) + ... + A\X'(t) + AoX(t) = 0 mi i (1) 

with known initial conditions 

X(t 0 ),X'(t 0 ),...,X n -\t 0 ) (2) 

where A,, i = 0,1,..., n G Af (to x r;X), (i.e. the algebra of square matrices with elements 
in the field X) with X(t) G M{m x 1; X) and detA n =0 if A n is square. For the sake of 
simplicity we set M m = Mirn x to; X) and A4 mr = M{rn x r; X). In the sequel we adopt 
the following notations 

Y 1 (t)=X(t), 

Y 2 (t)=X'(t), 

(3) 

Y n (t) = X^-'Xt). 
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(4) 


Y{(t) = X'(t) = Y 2 (t), 

Y'(t)=X"(t)=Y 3 (t), 

Y n -. 1 (t)=X( n ~V(t) = Y n (t), 

A n Yn(t) = A n xW(t) = —A n -iY n (t) - ... - AiY 2 {t) - A 0 Fi(t). 

Or in Matrix form 

FY'(t) = GY(t) (5) 

with known initial conditions 

Y(t 0 ) (6) 

where Y(t) = [Y X T (t)Y^(t)... Y^ (t)] T (where () T is the transpose tensor) and the coeffi¬ 
cient matrices F, G are given by 



4n 

0m,m 

••• 0m,m 

0m,m 


0m,m 

• • • 

0m,m 


0m,m 

-M71 

••• 0?n,m 

0m,m 


0m,m 

0m,m • • • 

0m,m 

F = 

0m,m 

0m,m 

••• Im 

0m,m 

G = 

0m,m 

0m,m • • • 

I'm 


0m,r 

0m,r 

••• 0m,r 

An 


o 

A\ ... 

A n —i 


with corresponding dimension of F, G and Y(t) mn x (mn + r — m), mn x (mn + r — in) 
and (mn + r — m) x 1, respectively. The matrix F is singular if A n (and consequently F) 
is singular. 

The matrix pencil theory has been extensively used for the study of Linear Matrix 
Differential/Difference Equations with time invariant coefficients, see for instance [1-6, 8, 
10, 11, 13-16, 19-23, 27-29, 36-42]. Recently this has been extended also to Linear Matrix 
Fractional Differential/Difference Equations, see [21, 22, 30, 31]. A matrix pencil is a 
family of matrices sF-G, parametrized by a complex number s, see [5, 6, 9, 20]. When G 
is square and F = where I n is the identity matrix, the zeros of the function det(sF-G) 
are the eigenvalues of G. Consequently, the problem of finding the nontrivial solutions of 
the equation 

sFX = GX (8) 

is called the generalized eigenvalue problem. Although the generalized eigenvalue prob¬ 
lem looks like a simple generalization of the usual eigenvalue problem, it exhibits some 
important differences. 

In the first place, if F, G are square matrices (for r=m) it is possible for F to be singular, 
in which case the problem has infinite eigenvalues. To see this, write the generalized 
eigenvalue problem in the reciprocal form 

FX = s~ 1 GX (9) 

If F is singular with a null vector X, then GX = 0 m7l; i, so that X is an eigenvector of 
the reciprocal problem corresponding to eigenvalue s -1 = 0; i.e., s = oo. In this case the 
solutions of system (1) have been fully analyzed in [1-3, 5, 6, 14, 19, 20, 22, 23]. 

Second, it is possible for det(sF-G) to be identically zero, independent of s. Finally it 
is possible both matrices F, G to be non square (for r^m). In this article we will consider 
these last two cases. 
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2 Singular matrix pencils: Mathematical back¬ 
ground and notation 

In this section we will give the mathematical background and the notation that is used 
throughout the paper 

Definition 2.1 Given F,G £ A i nm and an indeterminate s £ F, the matrix pencil 
sF — G is called regular when m = n and det (sF — G) ^ 0. In any other case, the pencil 
will be called singular. 

Definition 2.2 The pencil sF — G is said to be strictly equivalent to the pencil sF — G 
if and only if there exist nonsingular P £ Ml n and Q £ M m such as 

P(sF - G)Q = sF -G. 

In this article, we consider the case that the pencil is singular. The main results for this 
case are analytically presented. Unlike the case of the regular pencils, the characterisation 
of a singular matrix pencil, apart from the set of the determinantal divisors requires the 
definition of additional sets of invariants, the minimal indices. The distinguishing feature 
of a singular pencil sF-G is that either mn ^ mn or mn=mn and sF-G= 0. Let M r , Mi 
be right, left null space of a matrix respectively. Then the equations 

(sF - G)U(s) = 0 m ,i (10) 

V t (s)(sF-G) = 0 1iTO (11) 

have solutions in U(s), V(s), which are vectors in the rational vector spaces M r (sF — G) 
and Mi(sF — G) respectively. The binary vectors X(s) and Y T (s) express dependence 
relationships among the colums or rows of sF-G respectively. U(s), V(s) is polynomial 
vectors. Let d=dimA/ r r (sF — G ) and t =Mi(sF — G). It is known [5, 6, 9] that M r {sF — G), 
Mi{sF—G), as rational vector spaces, are spanned by minimal polynomial bases of minimal 
degrees 
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respectively. The set of minimal indices e* and Q are known [5, 6, 9] as column minimal 
indices (c.m.i.) and row minimal indices (r.m.i) of sF-G respectively. To sum up in the 
case of a singular matrix pencil, we have invariants, a set of elementary divisors (e.d.) and 
minimal indices, of the following type: 

• e.d. of the type (s — a ) Pj , finite elementary divisors (nz. f.e.d.) 

• e.d. of the type s q = -Y, infinite elementary divisors (i.e.d.). 

• m.c.i. of the type ei = 62 = ... = e g = 0 < e g +i < ... < Cd , minimal column indices 

• m.r.i. of the type = C ,2 = ••• = (h = 0 < Ch+i < ••• < Ct> minimal row indices 
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Definition 2.3. Let Bi, B 2 ,..., B n be elements of M n . The direct sum of them denoted 
by Bi ® B 2 © ... © B n is the blockdiag [ Bi B 2 ... B n ]. 


The existence of a complete set of invariants for singular pencils implies the existence 
of canonical form, known as Kronecker canonical form [5, 6, 9] defined by 

sFk — Qk sip — J p © sH q — I q ® sF e — G e © sF^ — G f © 0h, g (14) 

where sl p — J p , sH q — I q are defined in section 2. Let B 1 , B 2l ..., B n be elements of M n . 
The direct sum of them denoted by Bi © B 2 © ... © B n is the block diagjBi, B 2 , ..., B n }. 
where sl p — J p is uniquely defined by the set of f.e.d. 

V 

(s-ai ) P1 ,..., (s - a v ) p ", ^2 Pj =p (15) 

i =1 


of sF-G and has the form 


sip J p .— slp^ Jpi (^ 1 ) ® * * * © slp u Jpv (^^) (1^0 

The q blocks of the second uniquely defined block sH q — I q correspond to the i.e.d. 

(7 

si\...,s q °, £<?,' = 9 ( 17 ) 

i =1 


of sF-G and has the form 


sH q Iq :— sBq 1 I qi © ... © sHq a Iq a (IS) 

Thus, H q is a nilpotent element of M n with index q = max{<jj : j = 1,2,..., cr}, where 

H'l — f) 

and I p ., J Pj (a,j ), H q . are defined as 

'1 0 ... 0 O' 

0 1 ... 0 0 

Ipj — ; . . . € A4 Pj , (19) 

_ 0 0 ... 0 1 _ 

a,j 1 ... 0 0 

0 cij ... 0 0 

Jpj( a j) = : : : : & Mpj (20) 

0 0 ... aj 1 

0 0 ... 0 a :j _ 

'0 1 ... 0 O' 

0 0 ... 0 0 

- ; ; ■•. ; ; eM qj . ( 21 ) 

0 0 ... 0 1 

_ 0 0 ... 0 0 _ 
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For algorithms about the computations of the jordan matrices see [5, 6, 8, 9, 11, 18, 24, 
25]. For the rest of the diagonal blocks of Fk and Gk , s F e — G e and s F^ — Gq , the matrices 
F e , G e are defined as 




F e = blockdiag {L eg+ 1 ,L €g+2 ,. 


(22) 

Where L e = 

. h 

: 0 £>1 

, for e = e g+ 1 ,...,e d 





G e = blockdiag {L eg+1 , L £g+2 ,. 


(23) 


Where L f = 


Where Lq = 


Where Lq = 


1 '■ / e J j f° r e = € g+ 1 ) •••> The matrices FG q are defined as 

F ( = blockdiag {L Ch+1 , L (h+2 ,..., L Ct } (24) 

, for C = Ch+i, • • • > Ct 

G c = blockdiag {1^,1^, ...,L Ct } (25) 

°r' < ’ , for C = Ch+l)->Ct 
k . 


k 

Oi,C 


3 Solution space form of LMDEs with singu¬ 
lar matrix pencil 

Following the above given analysis, there exist non-singular matrices P, Q such that 

PFQ = F k 



PGQ = G k 

(26) 

Let 

Q — [ Qp Qq Qe Qc, Qg ] 

(27) 

where Qp (E Qq 

^ ^(mn)^ Qe £ Q £ ^ • / ^d(mn)£ ^^d Qg £ mn)g 


Lemma 3.1. System (5) 

is divided into five subsystems: 



io 1 

N? 

II 

HO 

(28) 

the subsystem 

H q Z' q {t) = Z q (t) 

(29) 

the subsystem 

F e Z' e (t) = G e Z e (t) 

(30) 

the subsystem 

F ( Z' ( (t) = GMt) 

(31) 
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the subsystem 


( 32 ) 


0 h,g ' Z g {t) — 0 h,g ' Z g (t) 

Proof. Consider the transformation 

Y(t) = QZ(t ) (33) 

Substituting the previous expression into (5) we obtain 

FQZ'{t) = GQZ(t). 

Whereby, multiplying by P and using (14), (26), we arrive at 

F K Z'{t) = G K Z(t). (34) 

Moreover, we can write Z(t) as 

Z p (t) 

z q {t) 

Z(t) = Z e (t) 

Z&) 

. z g (t) _ 

where Z.p^t') G .Adpi, Zq(t ) G •Mqi, Z e {t^j G .Adpi, Z^(t.') G J^dq\ and Z g (t ) G .Adpi, 
Z q {t) G M q \. Taking into account the above expressions, we arrive easily at (28-32). 

Solving the system (5) is equivalent solving subsystems (28-32). 

Remark 3.1. System (28) is a regular type system and its solution is given from, see 
[1-3, 7, 8, 10, 19, 20] 

Z p (t) = e Jp(*-*o)z p (t 0 ) (35) 

Remark 3.2. System (29) is a singular type system but its solution is very ease to 
compute, see [5, 6, 14-16, 22, 23] 

Zq(t)=0 qA (36) 

Proposition 3.1. The subsystem (30) has infinite solutions 


Proof. If we set 

Ze(t) 

The system (30) can be written as: 


Ze g+1 (t) 

Z eg+2 {t) 

Ze d {t) 


T Z 'e g+ ,(t) 
Z' (i) 

blockdiag{L eg+1 ,...,L Cd } 

- Z 'e d (t) 


= blockdiag {L eg+1 , 


Z £g+i ft) 

j Ze e+2 (t) 
Ze d (t) 
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Then for the non-zero blocks we have: 


L e Z' e (t) — L e Z e (t) , e — e g+ i,ed 
Using (22), (23) a typical equation from (37) can be written as 


h 


0 


e.l 


m = 


0 


e,l 


Ze{t) 


(37) 


(38) 


Every solution of the system (37) can be assigned arbitrary. The system (38) is of a regular 
type differential system. It is clear from the above analysis that in every one of the d-g 
subsystems one of the coordinates of the solution has to be arbitrary by assigned total. 
The system has no unique solution and can be taken arbitrary 

Z e (t) = C kA (39) 

Proposition 3.2. The subsystem (31) has the unique solution 

Zdt)=%A (40) 


Proof. If we set 



Z c {t) = 

Zch+ iW 
Zch+ 2W 



The system (31) can be written as: 

. z Ct (t) _ 



blockdiag {L Ch+1 ,L (t } 

i 

+ + * * 

N N 

1 _ 

= blockdiag {L (h+1 ,L^} 

’ z Ch+ 1 {t) ' 

z ( h +2 (*) 


L z iw \ 



. z Ct (t) _ 

Then for the non-zero blocks we have: 





L C Z' c {t) - L c Z c {t) , C — Ch+i, Ct (41) 

Because of structure of (L^,L^) blocks, it is readily shown that the only solution of (41) 
is the zero solution. 

Z c (t)=0 t -h,i (42) 

Remark 3.3. The last system (32) has an infinite number of solutions that can be taken 
arbitrary 

Z g (t) = C k , 2 (43) 

Theorem 3.1. Consider the system (5), with known initial conditions (6) and let the 
matrix pencil sF-G be singular. Then the solution is unique if and only if the c.rn.i. are 
zero 

dimN r (sF — G) = 0 (44) 
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and 


(45) 


Y(to) €E colspanQ p 

The unique solution is then given from the formula 

Y(t) = Q p e Jp( - t ~ t °' > Z p (t 0 ) (46) 

where Z p (to) is the unique solution of the algebraic system Y(to) = Q p Z p (to)- In any 
other case the system has infinite solutions. 

Proof. First we consider that the system has non zero c.rn.i and non zero r.nr.i. Consider 
the transformation (33), then from lemma 3.1, remarks 3.1, 3.2, 3.3 and propositions 3.1, 
3.2, the system (5) is divided into the subsystems (28-32) with solutions (35), (36), (39), 
(40), (43) respectively, then 


Y(t) = QZ(t) = [Q P Q q Q e Q ( Q g 


e Jp(t to)Z p (to) 

0,.i 

Ck, 1 

0t-h,l 

Ck, 2 


Y(t) — Q p e Jp ^ t to ^Z p (to) +Q e Ck,i + Q g C ka 

Since C k ,i and C k ,2 can be taken arbitrary, it is clear that the general singular LMDE for 
every suitable defined initial condition has an infinite number of solutions, so it doesn’t 
represent a dynamical system. It is clear that the existence of c.rn.i. is the reason that tlis 
system (30) and consequently system (32) exist. These systems as shown in propositions 
3.1 and remark 3.3 have always infinite solutions. Thus a necessary condition for the 
system to have unique solution is not to have any c.rn.i. which is equal to 

dimJ\f r (sF — G) = 0 

In this case the Kronecker canonical form of the pencil sF-G has the following form 

sFk — Qk := sip — J p ® sH q — I q ® sF ( * — G £ (47) 

and then the system (5) is divided into the three subsystems (28), (29), (31) with solutions 
(35), (36), (40) respectively. Thus 


Y(t) = QZ(t) = [ Q P Q q Q c 


e Jp(t *o)z p (t 0 ) 

0q,l 

0 t-h,l 


Y(t ) = Qpe Jp ^-^ Z p (to) 

The solution that exists if and only if 


Y(t 0 ) = Q P Z p {t 0 ) 



or 


Y(to) £ colspanQ p 

In this case the system has the unique solution 

Y(t) = Q p e J ^- to) Z p (t 0 ) 


Remark 2.2.5. It follows that the system (1) with known initial conditions (2) has a 
unique solution if and only if the initial value problem (5), (6) has a unique solution. Then 
its analytic solution is given by 

X(t) = Qy^z p (to) (48) 

where Q * is defined from the matrix 


Qp 


Ql ' 
Ql V 


Qp £ Mpp and where Q 1 is defined from the matrix 


and Q £ ^f(mn)(mn) 



Q 1 ' 
Q 2 


4 Numerical example 

4.1 Example 1 


Consider the matrix differential equation (5) and let 
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0 
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0 
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0 

0 

0 

0 
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0 

0 

0 

0 

0 
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Then det[sF-G]=0. The invariants of the pencil are, s-2, s-1 the finite elementary divisors, 
£i=0, €2=2 the c.m.i and Ci=0i C 2 =l are the r.rn.i. From theorem 3.1 the solutions of the 
system are infinite. 
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4.2 Example 2 

Consider the matrix differential equation (5) and let 

'11111' 

0 110 1 
11111 
0 110 1’ 

10 10 0 

0 0 111 

and 

12 2 12 

0 2 2 0 2 

1 2 2 2 3 

0 2 3 1 3 

0 0 0 0 0 

10 10 0 

The matrices F, G are non square. Thus the matrix pencil sF-G is singular with invariants, 
s-2, s-1 the finite elementary divisors, the infinite elementary divisors are of degree 1 and 
Ci=0, £ 2=1 are the r.m.i. From theorem 3.1 there exist non singular matrices 

'1 -10000' 

0 1 0 0 0 0 

-1 0 1000 
0 00010 
0 0 0 0 0 1 

0 -10100 

and 

'001 1-1 
1 1-1-1 0 

Q = 0 0-10 1 

1 0-1-1 1 

-10 2 1-1 

such that 

1 0 0 0 0 

0 10 0 0 

0 0 0 0 0 

0 0 0 1 0 

0 0 0 0 1 

0 0 0 0 0 

1 0 0 0 0 ] 

0 2 0 0 0 

0 0 10 0 

0 0 0 0 0 

0 0 0 1 0 

0 0 0 0 1 



PGQ = G k = 


PFQ = F k = 


10 



with 


0 0 
1 1 
Q P = 0 0 

1 0 
-1 0 

Let the initial values of the system be 

0 

-1 

y(o) = o 

i 

-i 

Then 

F(0) £ colspanQp 

and from theorem 3.1 the solution of the system is 

0 0 
1 1 
0 0 
1 0 
-1 0 

and by calculating Z p ( 0) we get 

'0 0 " 

1 1 

F(0) =00 Z p ( 0) 

1 0 
-1 0 

or 

z P { o) = _ X 2 

and the solution of the system is 



Next assume the initial conditions 

' 0 ' 

0 

Fo = 0 

1 
1 




(50) 
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Then 


To ^ colspanQp 


the initial conditions are non consistent and the solution for every k > 0 is 


Y{t) 


0 0 

e t e 2t 

0 0 

e t 0 

-e* 0 


C 


where C= 
is 2. 


Cl 

C2 


is constant and the dimension of the solution vector space of the system 


Conclusions 

In this article, we study an initial value problem of a class of linear rectangular matrix 
differential equations. First by taking into consideration that the relevant pencil is singular, 
we decompose the autonomous linear differential system into five sub-systems and we 
provide necessary and sufficient conditions for existence and uniqueness of solutions. As a 
further extension of the present paper, we can discuss the non-homogeneous case and the 
case of a singular discrete time system with a singular matrix pencil. For all this there is 
some research in progress. 
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